ISOMETRIC EMBEDDINGS IN IMAGING AND VISION: 
FACTS AND FICTION 



EMIL SAUCAN 

Abstract. We explore the practicability of Nash's Embedding Theo- 
rem in vision and imaging sciences. In particular, we investigate the 
relevance of a result of Burago and Zalgaller regarding the existence of 
isometric embeddings of polyhedral surfaces in and we show that 
their proof does not extended directly to higher dimensions. 



1. Introduction 

Recently, the attention of the vision community was drawn towards the 
important problem of isometric embeddings of manifolds in ([43]), in 
particular to the classical case of surfaces and their embedding in or R^ 
([9]). This approach was further extended by the employment of metric tech- 
niques that achieved fame through the celebrated work of Michael Gromov 
[2U] - see, [9], [10], amongst others. We do not consider here the critique 
of the metric method in general, since that would be, indeed, against our 
belief that it represents a basic, essential tool in imaging (see, for reference, 
[62]), but concentrate rather on the first problem: that of the isometric 
embedding. 

Let us begin by noting that, certainly the departure from the entrenched 
view of surfaces as given only locally by parameterizations (usually via by 
spline functions) has merit, as well as the realization that surfaces, that is 2- 
dimensional manifolds need often to be embedded, in a metrically controlled 
manner, m some R^, the value of N being, preferably, equal to 3. However, 
there are some deep and disturbing problems stemming from this approach, 
that is based upon the celebrated Nash Embedding Theorem. 

We begin with the least of these problems: the feasibility of finding an 
isometric embedding of a given, smooth orientable surface (or, more gener- 
ally, an orientable manifold) in some R^, for large enough. The root of 
the difficulty in writing and implementing an algorithm based upon Nash's 
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Theorem resides in the fact that this theoreni^ is a obtained via a fixed 
point method. The impediment here resides not only in the fact that such 
a method requires, at least theoretically, an infinite number of step^, but 
rather in the disturbing fact that the manifolds in the approximating se- 
quence are not usually submanifolds of the same M'^, where N represents 
the dimension of the Euclidean space in which the target (approximated) 
manifold, is (ideally) to be embedded - see also Remark 1.1 belowH This is 
not just a theoretical, quasi-philosophical quandary, but rather a serious im- 
pediment: indeed, since computers, by their very nature, can perform only a 
finite number of approximations, the computed manifold will be embedded 
in a dimension different from the one stipulated (and often depicted - see 
also Remark 15.161 below) . 

Even if we can, somehow, surmount this difficulty, one still is faced with 
the dire specter of dimensionality. To make this assertion clearer (and more 
concrete) let us recall a few facts regarding Nash's Theorem: 

The celebrated Nash embedding theorem assures the existence of an iso- 
metric embedding of any C^, {3 < k < oo) orientabl^ manifold of dimension 
n into some M^, for some sufficiently large. In contrast with the topo- 
logical case, were one has the classical Whitney's Theorem (see Appendix 
1), the embedding dimension ensured by the Nash theorem is prohibitively 
large, giving N = Zii^lLhlll^ (Note, however, that in Nash's original work, the 
embedding dimension for noncompact manifolds was N = |(n2+5n)(n+l).) 

Even with the further dimension reductions of Gromov ([26j. \28\ ) and 
Giinter ([30], EI]), the assured embedding dimension, namely (n^ + 10n-|-3), 
for 71 > 3, and (n + 2)(n + 3), if n > 4 and (max { "^"^^^ + 5, "("+^) |^ j-g. 
spectively, numerics are not very promising: the embedding dimension of 
a surface (i.e. 2-dimensional manifold) provided by the original Nash The- 
orem is 17, and by Gromov's and Giinter's improvements, 10. However, a 
special method developed by Gromov (|28], p. 298) decreases the embed- 
ding dimension for compact surfaces to 5, while for compact 3-manifolds 
the lowest guaranteed embedding dimension is = 13 - see [28j, p. 305. 
(Moreover, there exists a local isometric embedding of a given in any 
5-dimensional manifold - see [28j.) 

Strikingly, there is nothing really known about the case k = 2 (its omission 
in the discussion above being no mistake) Since curvature (of differentiable 



together with other two recently famous methods: the circle packing ([39], [40]) and 
the Riemann Mapping Theorem (see, e.g. ^). 

^ and that it is far less algorithmic in nature than the Picard Fixed Point Theorem of 
Differential Equations fame 

very similar problem arises also in the computational Riemann Mapping Theorem. 
^In the following, all manifolds are supposed to be orientable, except if otherwise specif- 
ically stated. 

^it is not even known whether manifolds admit isometric immersions in R'^ 
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surfaces, at least) is essentially a notiorH one would count on a general 
Nash embedding theorem for this case, that would allow a straightforward 
application in imaging. 

When one further relaxes the smoothness condition, the embedding di- 
mension decreases dramatically: any orientable 2-manifold is isometri- 
cally embeddable in M.^ (see ^0\, [55] )Q Indeed, since the proper (i.e. dif- 
ferentiable) notion of curvature makes sense only for manifolds of class > 2, 
the additional dimensions that are required to deal with curvatureH are not 
necessary in the case. The role of the lack of differentiability is further 
emphasized by the fact that by Nash-Kuiper Theorem, even the flat torus 
is isometrically embeddable in W^, in contrast with Tompkins' Theorem 
[70] (see also [65], pp. 196-197) that asserts that compact flat manifolds 
are not even isometrically immersable into M^"^^. Not only this, but also 
the following result (due to Kuiper [38]) holds: The unit sphere S" C W^~^^ 
admits an isometric immersion in M""*"^, for any n > 1. (This being in 
sharp contrast with the fact that any such immersion is congruent to the 
unit sphere, for n>2.) 

Moreover, any embedding is smoothable to C^, therefore Whitney's 
Embedding Theorem (see Appendix 1) implies Nash's Embedding Theorem 
(for n > 2). 

We should note that the case of analytic manifolds is, again, different: In 
[52j Nash proved that any compact (real) analytic n-dimensional manifolds 
has an isometric embedding in |(n^ -|- 5n). Gromov [26] extended Nash's 
result to include noncompact manifolds and also reduced the embedding 
dimension (for both cases) to i(n2 -h7n-M0). 

Remark 1.1. The striking difference between the results for the various de- 
grees of smoothness emphasize once more the delicate manner in which one 
should approach the embedding problem, and even more so its practical ap- 
plications. Moreover, even small variations of the metric, especially those 
producing change in the sign of (Gaussian) curvature (see also Remark l5.2p . 
can abruptly chang^ the embedding dimension (see [33] for a plethora of 
results in this direction). 

It follows that, when using approximating sequences of PL manifolds 
and the isometric embedding technique, one can not ascertain with any 
degree of certainty that the sequence of embeddings remains in the same 
(minimal) dimension. Indeed, those familiar with polygonal meshes (e.g. 
people working in the field of graphics) know that - quite counterintuitively 



Indeed, almost all the corpus of classical differential geometry of surfaces may be 
developed assuming only this degree of smoothness. 

'''The PL version of this result (and its extension to higher dimensions) represents the 
subject of following 3 sections. 

^using, for instance, the Gauss Equation - see, e.g. [T], [18] 

^and, unfortunately, usually increase 
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- even polygonal approximations of spheres have vertices of (concentrated) 
negative curvature (that is saddle points). 

Therefore, the only general assured embedding dimension is that guaran- 
teed by the Nash-Gromov-Giinther Theorem. 

Remark 1.2. Since we perceive shape, rather than distance, and since, by the 
previous remark (and the discussion preceding it, as well as by Remark 3.5 
below), curvature, hence shape, is lost, it follows that isometric embed- 
dings and their PL counterparts are far less useful in imaging, recognition 
and matching purposes than hoped for. 

Before proceeding any further, we have to add that we are aware that 
some readers are less familiar with some of the necessary background in 
differential topology and that, in any case, it would be best to refresh the 
basic necessary notions regarding immersions and embeddings. Therefore, 
we have included a glossary of relevant notions as Appendix 1. Given its 
goal, the pace is slow and the tone is rather didactical, therefore many a 
reader may want to omit it. However, since we also discuss therein one of 
the common misunderstandings regarding isometric embeddings, we believe 
it may be useful to all. 

2. A SUGGESTED SOLUTION: PL ISOMETRIC EMBEDDINGS 

It has been suggested to us [57] that the disquieting facts regarding the 
smooth embeddings considered in Nash's Theorem need not disturb us to 
much, for, in the imaging and graphics practice, one is always faced with 
PL-flat surfaces ("triangular meshes") and for these a highly surprising 
and widely unexpected result exist^, namely the following theorem due to 
Bur ago and Zalgaller: 

Theorem 2.1. Any compact orientable PL 2-manifold admits an isometric 
embedding in M^. 

The common wisdom regarding the statement above is, of course, that 
in imaging and graphics such surfaces represent the geometric object under 
investigation, or at least a "decent" approximation of it. 

Unfortunately, this does not represent the solution of the problem in ques- 
tion. Indeed, a number of problems arise as soon as one examines this 
theorem a bit closer. 

To begin with, the formulation above, while convenient and easy to recall, 
is not the correct one. The correct one can be found in [28j, p. 213 (but 
recall also the title of the paper [15j of Burago and Zalgaller): 

Theorem 2.2 (Burago-Zalgaller [28]). Every compact oriented surfaces with 
a piecewise linear metric can be piecewise linearly isometrically embedded in 

M3. 



'For some more recent, seemingly paradoxical, related results, see e.g. [13], [53], [54| . 
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First thing that strikes us is the apparently cumbersome and futile new 
terminology. However - as usually is the case in Mathematics, these apparent 
pointless minutiae and stresses are essential, and not due to just a whim 
of the mathematician. To comprehend this better in the case at hand, we 
should first understand the difference between the two (apparently identical) 
notions 

PL isometric embedding If P is a simplicial polyhedron of dimension 
n, any simplicial map f : P ^ (which is, by definition, linear on any 
simplex of P) induces a flat metric on each such simplex and, in consequence, 
a (singular) Riemannian metric g on P. More precisely, if P has k edges, 
then g is uniquely determined by the the vector {gi, ■ ■ ■ , gu), where gi = 
(length(ej))^. Saying that / is a PL isometric embedding means that / is 
as above and, in addition, it is also an embedding. 

PL isometric embedding of subdivided polyhedra In this case, the map- 
ping (embedding) f : P ^ M"* is required to be linearly isometric on the 
simplices of a simplicial subdivision of P. Evidently, by a sufficient num- 
ber (albeit practically infinite) number of subdivisions, one can approximate 
the Riemannian case using PL metrics, i.e. such that each simplex of P is 
isometric to a Euclidean simplex (in R"). 

So, why is the first definition not adequate? The main problem is its 
rigidity: Informally put, one "has to work with what he's got". That is, 
further subdivisions (hence approximations) are not allowed. Therefore, 
this approach rapidly reduces to a largely combinatorial problem, at least 
in many of its aspects (see, e.g. [37]). 

Moreover (and more important) this rigidity is not just of convenience 
(so to say), quite the contrary - it is essential. Indeed, mos PL isometric 
maps are rigid, in the (geo-)metric sense: 

Theorem 2.3 ([28j). Every small deformation of a n-dimensional polyhe- 
dron embedded (immersed) in R"+^ is an isometry. 

In fact, the result just mentioned is more general, but to avoid a further 
detour, we refer to [28]. pp. 210-211.) 

Obviously, the second notion is far more attractive, both for the geome- 
ter/analyst as well as for Computer Graphics and related fields. However, 
caution should be taken, since this is still a very fiexible notion, since just 
the metric is to be preserved. For instance, one has the following result of 
Zalgaller ([73]): 



Given the space limitations and the desire for cohesiveness, we must assume the 
reader is familiar with the very basic notions of PL topology. (For a deep, yet enjoy- 
able and not overly technical source on these notions, see [69] . See also Appendix 1 for 
additional material on embeddings.) 

Here "most" has a precise mathematical meaning: More precisely, generic simplicial 
mappings are rigid, where "generic" is a rather technical term (see [28]). 
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Theorem 2.4 (Zalgaller, 1958). Let P he a simplicial polyhedron of dimen- 
sion n,n < 4, endowed with a PL metric. Then P admits an equidimen- 
sional PL map into M". 

This is a very surprising and counterintuitive result 1^ So, even though 
the embedding condition is omitted, it prepares us to understand somewhat 
better the problematic nature of the definition and of the Bur ago- Zalgaller 
theorem. 

To further elucidate these notions we take below a closer look at Burago 
and Zalgaller's proof. 

3. The Burago-Zalgaller gonstruction 

Main idea of the proof is - not very surprisingly - to adapt the proof 
of the Nash- Kuiper Embedding Theorem^ For this, one starts with 
a (smooth) embedding of the given polyhedron in M'^, composed with a 
contracting homotety. Then: 

• Carry out a sequence of stages, divided in turn into a large number 
of steps, each of which improves the approximation to isometry and 
such that the function obtained at each stage is short: 

Definition 3.1. {Short mappings) Let {X^d) and iY,p) be metric 
spaces. A map / : X — )• y is called C-short iff 

p{f{x), f{y)) < Cd{x, y) , for allx,yeX. 
f is called short (or contracting) iff it is C-short, for some C < 1. 

• Add "ripples" , producing thus a PL version of Kuiper's adaptation 
|38j of Nash's twist [50j) such that one will have "enough space" to 
isometrically embedded the surface in M'^. 

We won't dwell too much in the details of the proof, just mention some 
of the principal "geometric" stages: 
(1) Basic Construction Element 

(a) Let T = A(^i, A2, A3) and t = A(ai, 02, 03) be acute triangles; 

(b) let B,b and R,r the centers and radii of their respective cir- 
cumscribed circles; 

(c) let Ep = \AkAi,ep = ^auai; p,k,l £ {1,2,3}; 

(d) and let Hp = BEp, hp = bcp. 

Moreover, let T ~ AkAi > a^a/, k, I G {1, 2, 3}. 

Then T can be isometrically PL embedded in M^, as the pleated 

surface included in the right prism with base t, such that AkAiAp 

^"^However, it represents the PL version of the version (that holds for any n) - see, 
e.g. [28]. 

^^We preserve in the following overview of the proof the notation of [15]. For copyright 
reasons we do not reproduce, however, the figures included therein, but rather refer the 
reader to the original source. 
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fits akEpttiE'^apE^, where: B'bJ-t, B'ap = R and Ep^E'^,E[ on 
the faces of the prism, such that aj^E'p = Eipai = ^Af^Ai. 
The following variations of the basic construction above are also 
considered: 

(i) Each angle (p oiT satisfies the condition {) < a < (p and 
C ■ AkAi > akCLi , C < 1. Moreover, AkAi/akai ~ 1. 

(ii) Each of the lateral faces of the prism - including the bro- 
ken lines akEptti - can be (independently) slightly rotated 
around the lines a/jO/ such that the construction still can 
be performed. (The rotation angle depends upon the con- 
stants a and C above.) 

(In general, one has to simultaneously construct a large number 
of the units above.) 

(2) Standard embedding near vertices 

Use the standard conformal map (or folding) from K{9, p) = {0 < 
ip < 9, p> 0} to K{X, r) = {0 < V < A, r > 0} given by: 

Tp = ^(p, r = ap^/^. 

(The most important case for our purposes being: A = 2ii.) 

(3) The Triangulation and its refinement 

Let ^4 be a vertex of total angle 9. 

(a) If 9 < 2tt, then encircle ^ by a small "regular" hexagon com- 
posed of 6 triangles of apex angle 9/6. 

Some small enough neighbourhood of A the will be mapped by 
the standard conformal mapping onto a planar disk. 
Over each triangle included in such a neighbourhood, one can 
perform the basic construction, obtaining a PL isometric em- 
bedding of this neighbourhood. 

(b) If 6* > 2tt, proceed analogously to the previous case but 

(i) In a small circular neighbourhood of radius ri map (a) 
isometrically on radial segments and (b) using a 9/2pi 
contraction on circles centered at A; 

(ii) In a annular neighbourhood {ri < r < use the stan- 
dard conformal mapping with the same contraction factor 
9/2pi. 

Replace the neigbourhood above with a "cogwheel" (i.e. a cir- 
cle surrounded by isosceles "triangles" of sides, e.g. 26, and 
having as bases arcs of the same length). The interior of each 
"cogwheel" is PL isometric embedded using "ripples". (The 
basic element of each such "ripple" is a pair of congruent trian- 
gles, having a common vertex in the center of the "cogwheel", 
one side (of each) being a radius, and a second common vertex 
built over the midpoint of an arc used in the construction of the 
"cogwheel" - see Figure 4 of plSj). Away from neighbourhoods 
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of vertices, refine the triangulation using only acute triangles. 
In particular, at convex vertices subdivide each triangle into 
similar triangles, for some large enough n; while at non-convex 
vertices into almost regular triangles. 



Remark 3.2. Adaptations [15] of the main technique exposed above ensure 
the existence of PL isometric embeddings of (orientable) PL manifolds with 
boundary and of PL immersions of nonorientable PL manifolds. 

Remark 3.3. A close examination of the arguments of the proof shows that 
Theorem 2.2 can be extended to include (orientable) non-compact manifolds 
with bounded (generalized) principal curvatures |60j . 

Yet one naturally has to ask himself the following 

Question 1. Is the Burago-Zalgaller Theorem applicable for Image Process- 
ing/Computer Graphics? 

Unfortunately, the answer is negative, for the following reasons: 

(1) The construction yields an (infinite) approximation process, akin to 
the original Nash-Kuiper method, hence numerical errors have to 
contended with and taken into account. 

(2) The geometrjEl of the limiting object is very far from the one of 
the "target surface": Not only is the resulting PL surface strongly 
"corrugated" (as evident from the construction), it my also contain 
"superfluous" vertices, i.e. where the curvature (of the metric) is 
zero. Moreover, for surfaces of positive extrinsic curvature, it is quite 
possible that the surface admits not even an isometric immersion in 

such that the extrinsic curvature equals the intrinsic one. (For 
the technical definitions and a simple example of a PL 2-sphere 
exhibiting this behavior, see [H], p. 76. 

Also, note that "accidents" in the original embedding can produce 
widely diverging subdivision schemes - see also the comment on page 
2 above, as well as Remark 1.1. (To grasp this widely divergent 
behavior, one should consider, for instance, the examples quoted in 
footnote 4, on page 2.) 

Example 3.4. {Burago-Zalgaller, Example 1.5) For any e > 0, the 
flat torus (i.e. the "topologist's torus", obtained by "gluing" the 
opposites sides of a plane square via Euclidean translations) admits 
a PL isometrical embedding e-close to the rotation ("round") torus 



all important in any practical implementation 

Again, this result represents the PL equivalent of its counterpart - see page 3 
above. 
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Remark 3.5. It is contended in [7J, p. 618, that the the Burago- 
Zalgaller embedding method_preserves curvature. However, this as- 
sertion is not made in [15jr1 Indeed, this is not possible, as the 
example above clearly hints and as we shall explain in some detail 
below. 

First, we should understand what type of curvature is preserved. 
Evidently, not the canonical ("smooth") one of classical differential 
geometry, since the considered surfaces are not even C^. It may be 
that the authors of [7] refer to the fact that a piecewise flat surface 
(or polyhedron) remains piecewise flat during the embedding process 
and, as such, its curvature is identically zero at all the points that do 
not belong to the vertices and edges of the triangulation. However, 
it is explicitly emphasized in the very introduction of |15j that "the 
metric of a polyhedron is locally flat excev r^ a t a finite collection of 
points; these points are the "true" vertices. 'o 

However, it is precisely at these points that Gaussian curvature is 
concentrated (being the defect of the planar angles (of the faces) in- 
cident at any such vertex). This is a known, in fact, since Descartes, 
but it was introduced in modern Mathematics by Hilbert and Cohn- 
Vossen [34], and developed first by Polya [58] and then by Banchoff 
[2], [3] and, more recently by Stone and Pu [21] 1^ (Mean cur- 
vature is, by contrast, concentrated along the edges of a polyhedral 
mesh, as the dihedral angle of the two faces who's intersection is 
any specific edge - see [l6] for a succinct presentation and for the 
bibliography within.) 



^^In fact, the authors of exphcitly state (see [H], p. 370) that they "... prove the ... 
discrete analog of the weU-known result of J. Nash and N. Kuiper on -smooth isometric 
immersions" (our emphasis), 
our emphasis 

^^[ig, p. 369 

^'^This definition of Gaussian curvature for polyhedral surfaces also facilitates an easy 
proof of the fact that no PL isometric embedding of the flat torus, preserving discrete 
curvature, is possible in R''. This proof follows closely the one for the smooth, classical 
case (see e.g. [IS]): Since the given polyhedral embedding V is compact, there exists a 
number R > such that V is included in the interior of §jj(0) - the sphere centered at the 
origin and of radius R. Let R decrease continuously and let _Ro be the radius for which the 
sphere and the polyhedron have the first non-void intersection. Moreover, by elementary 
arguments, any such intersection point must be a vertex of the polyhedron. Since V is 
contained in the closure of Sfjjj(O), that is on one side of the surface of the sphere, it 
follows that the sum of the angles at an intersection vertex vo must be smaller than 2n, 
that is the discrete Gaussian curvature K^g > 0, in contradiction to the supposition that 
P is a realization of the fiat torus, i.e. having discrete Gaussian curvature equal to zero 
at all its vertices. 



10 



EMIL SAUCAN 



4. A SHATTERED HOPE 

Having seen that the Burago-Zalgaller construction is not apphcable as 
such, one at least hopes for a positive answer to the fohowing natural 

Question 2. Does Burago-Zallgaler's Theorem hold in dimension n >3? 

Perhaps unexpectedly, and contrary to the unsubstantiated statement of 
00 the answer to this question is not known! However, there are indica- 
tions that the answer is negative. These indications emerge from the proof 
in dimension 2: 

(1) The proof is based on the previous result of Burago and Zallgaler on 
the existence of acute triangulations. 

Strangely enough, next to nothing is known about the existence 
of such triangulations in dimension n > 3p^ 

(2) The proof heavily relies on the use of the use of the standard con- 
formal map to produce a mapping that (around the vertices) is ar- 
bitrarily close to conformality (and, in the end) to isometry. 

However, in dimension n > 3 this is not possible: the analogue of 
the standard conformal map has dilatation bounded away from l! 
Indeed, we can be more specific. But first, a few technicao details: 

Definition 4.1. (Wedges) Let x G R" be a point with cylindrical 

coordinates x = (r cos (/?, r sin zi, . . . , z„_2). The set Da = {0 < 
ip < a}, (0 < a < 2tt) is called a wedge of angle a. 

Definition 4.2. (Foldings) The homeomorphism. / : Da — )• 
/(r. If, z) = (r, ^ip, z), z = (zi, . . . , Zn-2) is called a folding. 

Before proceeding further, the reader should familiarize herself/himself 
with the technical notions regarding quasiconformal mappings. Not 
wishing to interrupt the flow of geometric arguments, we have con- 
centrated these in a short appendix (Appendix 2). 

Proposition 4.3 (Gehring-Vaisala [2^,[7l]). Let Da, Di^ be wedges, 
and let f : Da — >■ Dp be the respective folding. If a < f3, then f is 
quasiconformal, with dilatations Kj(f) = ^,Ko(f) > (^)^^^"'~^^ . In 

particular, for (3 = vr, we obtain Ki(Da) = %Ko(Da) = (l)^'^'^~^\ 
whence K(Da) = ^■ 

Remark 4.4. Remarkably, the coefficients of quasiconformality for 
non-convex domains (i.e. vr < 2ti) are not known 



"^-'^In fact, the abstract ,44 does state the resuh but does not sustain it was actually 
obtained. 

*^^The little existing information is summarized by Zamfirescu [73]. (Some additional 

hope stems from a different method developed recently by Tasmuratov [67].) 
90 

But hopefully not too technical - we shall restrict ourselves to the simplest case. 

94- 

at least to the bet of our knowledge 
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Following |16] , we note the following natural generalization of the 
definition of a wedge: 

Definition 4.5. The domain D^k C W^,Dak = ■ ■ ■ ,^n-k+i, 

Zn-k+i, ■ ■ -jZn)}, < (pk < Ok ,1 < k < n-v-l,a = (ai, . . .,an~u~i), 
< ai < 27r, < 02, • • • , an-u-i < vr is called a dihedral wedge of 
type V and angle a. 

Remark 4.6. For k = n — 2 we recuperate the classical definition of 
wedges. 

The numbers that allow us to ascertain whether two domains are 
quasiconformally equivalent, i.e. that one is the the quasiconformal 
(therefore homeomorphic) image of the other, and, if so, which is the 
smallest possible dilatation of such a mapping, are called the coeffi- 
cients of quasiconformality (see Appendix 2). We have the following 

Proposition 4.7 ([16|). The coefficients of quasiconformality for 
Dak are: 



1 

^n— fc— 1 / ^n-fc-1 \ n-l 

Ki{Dak) = , Ko{Dak) > 



CHI ■ ■ ■ Un-k-l ■ ■ ■ an-k-l 

n-k-1 

K{Dak) 



"1 • • • a-n-k-i 

Corollary 4.8. Let V be a convex polyhedral domain in M" and let 
m denote the number of faces of V . Then we have the following 
estimates: 



m — n + 2 f m — n + 2\ "-1 

Ki{V) > , Ko{V) > 

m — n \ m — n 



m — n + 2 

K{V) > 



m — n 



Remark 4.9. Evidently, the same estimates hold for PL-smooth con- 
vex manifolds. 

Hence, for polyhedra with a very large number of faces, such as 
encountered in (good) PL approximations of domains D in (or, 
more generally, in M"), having smooth (convex) boundaries, K{D) ^ 
1. Even without considering approximations and without making 
appeal to Corollary 14. 81 one can easily produce (convex) polyhedra V 
that require arbitrarily large dilatation K{V), by choosing polyhedra 
with at least one dihedral angle (between n-faces) vr/m, where m is 
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any (arbitrarily large) natural number, and applying Proposition H31 
directlyEl 

It is interesting to note in this context, that (at least for polyhe- 
dral domains in M^) the "primary carrier" of dilatation is the mean 
curvature H - see Remark 13.51 above. 



5.1. Glimmers of hope. We bring below two different approaches to the 
embedding problem, that both circumvent the intricacies of the Nash Em- 
bedding theorem mentioned in the preceding sections. 

5.1.1. A compromise. Reviewing the facts above, it is hard not to reach the 
conclusion that the situation is quite bleak, as far as the practical use of 
Nash's Embedding Theorem is concerned. However, one may quite justifi- 
ably sustain that making appeal to global isometric embeddings in general, 
and to Nash's Theorem in particular, is to be somewhat overenthusiastic. 
Indeed, it may be very well claimed, that one is rarely faced, in computer vi- 
sion, graphics and other related domains, with surfaces (manifolds) globally 
defined, hence one can restrict himself to local isometric embeddings. After 
all, this is the position already adopted (albeit in a different context, where 
large amounts of data have to be processed) in the widely quoted work of 
Roweis and Saul [59] . The method of |68] is also basically local (see, however 
pO] for a discussion on its possible globality). 

This approach is also augmented by the very first result on isometric 
embeddings, namely the following theorem of Burstin, Janet and Cartan 
(see, e.g. [65]): 

Theorem 5.1 (Burstin- Janet-Cartan). Any (real) analytic manifold M" 

ii(n + l) 

can be locally (real) analytically isometrically embedded into M 2 

The problem with the result above is the fact that it requires analyticity. 
In fact, even if conjectured already by Schlaefiy in 1873, the proof of the 
result above for manifolds is still elusive. (It is true, however, that 
weaker forms of this result - that is, with higher embedding dimension - 



Remark 5.2. As in the global case, for lower differentiability classes, no 
general results are even possible. Indeed, there exist a counterexample, due 
to Pogorelov [56j, of a C^'^ metric on the unit disk = 1^(0,1) C M?, 
such that there exists no isometric imbedding in M.^ of B^(0,r), for any 
< r < 1. (See also [H], [l9] for some more recent results in this direction.) 



For a stronger result regarding the nonexistence of isometric embeddings for PL 
manifolds in dimension > 3, see [60] . 

The author is not aware of the existence of meaningful, general theorems regarding 
manifolds, for 2 < fc < 00. 



5. Final comments 
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Still, one may argue (rather convincingly) that it is quite commorP^ in 
imaging and vision to adopt smooth, even analytic models and consider stan- 
dard types of approximations (for the manifolds and for various differential 
operators on these manifolds). 

5.1.2. More dimensions. Surprisingly, a very effective (at least from the the- 
oretical viewpoint) alternative embedding method follows the quite opposite 
direction: Instead of reducing the scope of the embedding, one can extend 
it by adding dimensions. That is, one can embed not in some M'^, but 
in an infinitely dimensional space, more precisely in L°°(M") - the (Ba- 
nach) space of bounded Borel functions on M", endowed with the "sup" 
metric, i.e. d{f,g) = sup^^^" \f{x) - 9{x)\, for any f,g e L°°{M'^) - via 
the Kuratowski Embedding [35]: 

Definition 5.3. {Kuratowski embedding) Let be a closed Riemannian 
manifold. Then 

K -.W ^ L~(M") , K{x) = dist^ , 

where 

(5.1) dist^. = dist(x, •) , 

where "dist" denotes the (intrinsic, Riemannian) distance on M", is called 
the Kuratowski embedding (of M"). 

This method is much more powerful than it would appear at first sight. 
Indeed, the Kuratowski embedding is an isometry, more precisely we have 
the following Lemma (see, e.g. |32j): 

Lemma 5.4. With the notation above, we have 

d(dista;, disty) = dist(x,y) . 

Remark 5.5. This approach is widely divergent from the Riemennian embed- 
ding one adopted in Nash's Theorem. Indeed, the Riemannian and Kura- 
towski embeddings coincide iff K{M^) C L°^{M^) is a convex, open subset 
of an affine linear subspace of dimension n. 

On behalf of the Kuratowski embedding, one can remark that, albeit 
being infinite dimensional, it may be quite advantageous when a functional 
approach is needed or sought for (e.g. when considering spline functions, 
wavelets, etc.)J^ However, usually (and more realistically) the spaces that 
appear in Computer Science (and even more so in Graphics) are finitely 
dimensional. Moreover, most people in the said communities find infinitely 
dimensional spaces as somewhat of an artifice, highly nonintuitive, and of 
theoretical value at best. 

97 

even though the author does not subscribe himself to this philosophy 

closely related approach is well known to the Imaging and Vision community: Em- 
bedding by using the eigenvalues of the Laplacian or of the Green Kernel. (For applications 
of these methods in the context of Riemannian Geometry, see, e.g. [5], [41|.) 
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Fortunately, there exists a finitely dimensional version of the Kuratowski 
embedding: Let X be an e-ne^ in M", \X\ = m. Then, for small enough 
e, Kx : X ^ l'^ \s an embedding, where Kx = K\x - the restriction of K 
to X and denotes the m-dimensional Banach space endowed, again, with 
the "sup" metric: if x = (xi, . . . then ||x|| = supj \xi\. 

Moreover, we can assure that this "discrete" version of the Kuratowski 
embedding is bi-lipschitz, more precisely we have the following result ([32], 

my- 

Theorem 5.7. Let M" be a compact Riemannian manifold without bound- 
ary. Then, for any C > 0, there exists a e-net X, where e = e{C), such 
that 

(5.2) (1 - C)dist{x,y) < \Kx{x) - Kxiy)\ < dist(x,y). 

Due to the theorem above, the finite dimensional version of the Kura- 
towski embedding is proves to be very useful in Global Differential Geom- 
etry: Its use in the study of systoles was pioneered by Gromov [29j (see 
also [32], [l2]). It was also employed to prove yet another result of Gro- 
mov [27] (and Katsuda [37]), namely a rigidity theorem: Informally stated, 
the theorem in question asserts that if two re-dimensional (compact) Rie- 
mannian manifolds, having the same lower bound for their volumes, and 
upper bounds on diameters and sectional curvatures, are sufficiently close 
one to each other in the Gromov-Hausdorff (metric) topolog\P^. then they 
are diffeomorphic. 

Since the manifolds usually encountered in Imaging, Vision, etc., natu- 
rally satisfy such bounds, it follows that the result above, as well as the 
finitely-dimensional Kuratowski embedding in general, are quite relevant 
for applications in the mentioned fields, in particular for recognition type 
problems. 

5.2. A possible solution. A more realistic approach (both from the the- 
oretical and implementational viewpoints) would be to obtain a Discrete 
Nash Embedding Theorem [ITj. A certain amount of confidence in the fea- 
sibility of obtaining such a result stems, amongst others, from the existence 
of discrete versions of the required differential operators and invariants (see, 

e.g. m,m)- 

Another geometrization approach stems from the differential geometry 
of metric spaces (see, e.g. [8]). By using a discretization of the (metric) 
Finsler-Haantjes curvature of curves ([H]) we can obtain the embedding, via 
a proper discretization of the Gauss-Bonnet theorem, of any given metric 



Recall that e-nets are defined as follows: 

Definition 5.6. Let {X,d) be a metric space, and let A C X. A is called an e-net iff 
d{x,A) < e, for all x £ X. 

OQ |__| |__| 

See, e.g. [22], [62] for a short overview of the notion. 
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graph, not into M" (or Ip, H") but rather into a model space (a model surface, 
to be more precise) - see [63] . 

Naturally, one expects the two embedding methods considered above to 
converge and augment each other, particularly in our purely geometric con- 
text, stemming from problems in PL differential geometry, computer graph- 
ics and image processing, where the graphs considered are skeleta of trian- 
gulations of manifolds (or of cell-complexes), and the weights are either 
edge-weights (i.e. distances between vertices) or/and vertex-weights (i.e. 
curvature measures) - see, e.g. [62], [64] . 

However, both methods, applied in a more general context, rend them- 
selves to various practical implementations, in such areas as multicommod- 
ity flows in networks (e.g. for the prediction of informational bottlenecks, 
discovery of holes, etc.), clustering of statistical data (in particular in bio- 
informatics - see, e.g. [61]) and expanders. 
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Appendix 1 - Immersions and Embeddings 

We presume the reader is familiar with the notions of differentiable man- 
ifold and tangent space, as well as with basic concepts of topology (for any 
eventually needed details see ^]), and we recall only the relevant defini- 
tions: 

Definition 5.8. {Immersion) Let M^,N" be smooth differentiable man- 
ifolds and let / : M"^ N"' a differentiable map. If rank/ = m at each 
point of M"*, then / is called an immersion. 

Definition 5.9. (Embedding) Let M'^,N'^ be smooth differentiable mani- 
folds and let / : M*" — )• A^" a differentiable homeomorphism. If / is also an 
immersion, then it is called an embedding. 

(Note that, in this case, m = n, of course.) 

The condition that / be a homeomorphism is very strong and shouldn't 
considered lightly. Indeed, not even asking that / be injective will suffice, 
as proven by the (classical) fact that there exists a injective immersion of M 
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into the "figure eight" curve, but this is not an embedding, since it is not a 
homeomorphism: the image is not even a manifold 

To sum up: The notions of embbeding and immersion are not interchange- 
able - while any embedding is, in particular, an immersion, the opposite is 
not trueJB 

If one discards even the differential structure, then topological embeddings 
(in the sense of that they homeomorphic on their image) are relatively easily 
obtained by 

Theorem 5.10 (Whitney's Theorem). Every (smooth) manifold of dimen- 
sion n admits a (smooth) embbeding in M^" and a (smooth) immersion in 



Remark 5.11. It is easy to prove, for compact manifolds, that an embedding 
in some finite dimension exists. It is then progressively (much) harder to 
discard the compactness restriction and to gradually "zero in" to dimension 
2n, via embedding dimensions {n + 1)^ and 2n + 1. 

Up to this point we have dealt with classical ("pure") differential topol- 
ogy, that is the famous "rubber geometry" of popularization texts (albeit 
endowed with some "smoothness" - necessary for the "differential" part.) At 
this point, however, we should introduce a bit of "solid" geometry, neces- 
sary, e.g. for recognition purpose^. The idea is to use a specific "measuring 
yard" , for each manifold, that is a Riemannian metric: 

Definition 5.12. {Riemannian manifolds) Let M" be a manifold and let 
X = (xi, . . . , Xn) denote a standard coordinate chart. The Riemannian met- 
ric g on M" is defined by the length element 



i=l 

where the functions gij = gij{xi, . . . ,Xn) represent the scalar products of 
the vector fields -Mr associated to the given chart: 



Once an infinitesimal distance is introduced, global ones can also be mea- 
sured (transforming a Riemannian manifold into a "honest-to-God" metric 
space) as follows: 

Let c : [, ] — )■ M" be a curve. Then its length is given by: 

"'^An even more pathological example can be constructed, where R is injectively im- 
mersed in the (flat) torus as a dense geodesic (just consider the image under the 
covering map of T^, of a line making an irrational angle with the Ox axis. 

•^^Therefore, even though, for stylistic reasons, one usually tends to avoid repeti- 
tions of the same word, one cannot (see e.g. [10]) freely interchange "immersion" and 
"embedding" ... 

■^"^Think of the the individualization problems amoebas are faced with... 



n 
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length (c) = 




(Here is important to recall that c'{t) is just a tangent vector, so ||c'(f)|| = 
V5(c'(i),c'(t)).) 

The intrinsic (or inner) distance between two points p,q ^ M" is defined 

as 

d{p, q) = inf{length(c) | is a curve of ends p and q} . 

c 

(We have tried here to keep the technical aspects of the definition above 
to a minimal level; for those insisting on absolute formal correctness, we 
recommend, for instance, [17].) 

Of course, the Riemannian metric induces a topology on M" (the metric 
topology), but in fact a much stronger result holds: 

Theorem 5.13 (Palais, [55j). The metric of Riemannian manifolds deter- 
mines its (smooth) manifold structure. 

Thus, Riemannian manifolds are, a fortiori, smooth manifolds, and the 
discussion above holds for them as well. However, a much more specific 
notion of embedding is applicable (relevant) in this case, namely: 

Definition 5.14. {Isometric embedding) An embedding / : M" — >■ is 
called isometric iff 

(5.3) < >= gi,j A <hj <n, 
where < ■ , • > is the standard inner product in 

Remark 5.15. Of course one can extend the definition above to include em- 
beddings / : M" , where Q is a general Riemannian manifold with 
metric h, by imposing, instead of (|5.14p . the following condition: 

(5.4) < Vif, Vjf >h= QLj , 1 < i, J < n , 

where Vj/ = Df (^^^ ^''^^ where < •, • >h represents the scalar product 
defined by the Riemannian metric h on Tj(^^^{Q), (x = {xi, . . . ,x„) G M"). 

Remark 5.16. {A common fallacy) Sadly, the notion of isometric embedding 
(and in particular the result above) are sometimes puzzling even for the pro- 
fessional mathematician (and even, sometimes, for topologists!...) However, 
things become much simpler if one keeps in mind Definition l5.14l and remem- 
bers that saying that the intrinsic metric "equals" the Euclidean one, means 
just that the infinitesimal length element, as defined by the Riemannian 
metric coincides with that of the infinitesimal one induced by the ambient 
Euclidean space. In any case this shouldn't be interpreted as affirming that 
lengths of curves, as measured on the manifold, equal the Euclidean dis- 
tance between their ends, as measured in the ambient space (see Figure 
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below). Therefore, it would be redundant and useless (not to say mistaken) 
to try and isometrically embed such surfaces in , for some 




Figure 1. Intrinsic vs. ambient distance: The length of 
the (strait, added) cord never equals the one of the bow, 
independently of the dimension in which Rodin's bowman 
resides. 



On the positive side, surfaces in are already embedded and inherit, 
therefore, a Riemannian structure from the ambient space, that is the in- 
duced metric defines on the surface a Riemannian metric (even if none was 
supposed (or given) a priori). 



In fact, the only manifolds for which the embedding and intrinsic metric coincide are 
precisely the piecewise flat ones - see [6]. 
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6. Appendix 2 - Quasiconformal Mappings 

Definition 6.1. {Quasiregular and quasiconformal mappings) Let D C M" 
be a domain; n > 2 and let / : D — >■ M" be a continuous mapping. / is 
called 

(1) quasiregular (qr) iff 

(a) / is locally Lipschitz (and thus differentiable a.e.); 
and 

(b) < I/'(x)|" < KJf{x), for any x G M"; 

where = sup \f'{x)h\, and where Jf{x) = detf'{x); 

\h\ = 1 

(2) quasiconformal (qc) iS f : D ^ f{D) is a quasiregular homeomor- 
phism; 

The smallest number K that satisfies condition (b) above is called the outer 
dilatation of /. 

Remark 6.2. One can extend the definitions above to mappings between 
oriented, connected Riemannian n-manifolds, n > 2, by using coordinate 
charts (for details see, e.g. [7T|). 

Remark 6.3. It follows immediately from Condition (1) (b) above, that qr- 
mappings are sense preserving. 

If / : D — 7- M" is quasiregular, then there exists K' > 1 such that the 
following inequality holds a.e. in M": 

(6.1) Jf{x) < K' inf |r,.//ir 

\h\=l 

By analogy with the outer dilatation we have the following definition: 

Definition 6.4. [K -quasiregularity) The smallest number K' that satis- 
fies inequality (16. ip is called the inner dilation Ki{f) of /, and K{f) = 
max{Ko{f), Ki(f)) is called the maximal dilatation of /. If K{f) < oo we 
say that / is called K- quasiregular. 

The dilations are K{f), Ko{f) and Kj{f) are simultaneously finite or 
infinite. Indeed, the following inequalities hold: Ki{f) < K^^{f) and 
Ko{f)<K'}-\f). 

Definition 6.5. [Coefficients of quasiconformality) Let Di,D2 C M" be 

domains homeomorphic to each other. The numbers 

(6.2) 

Ko{Di,D2) = mfKo{f), Kj{Di,D2) = mfKi{f), K{Di, D2) = mf K{f), 

where the infima are taken over all the homeomorphisms / : Di — ?• D2 are 
called the outer, inner and total coefficient of quasiconformality of Di with 
respect to D2, respectively. If D2 is the unit ball B", then the numbers 
Ko{Di) = i^o(A,B"), KjiDi) = Ki{Di,W), K{Di) = K{Di,W) are 
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simply called the (inner, resp. outer, resp. total) coefficients of conformality 
of Di. 

Again, the numbers Ko{Di, D2), Kj{Di, D2) and K{Di, D2) are simul- 
taneously finite or infinite. However, it is not always guaranteed that there 
actually exists a homeomorphism / as above, such that Ki{f) = Ki{Di, D2) 
or Koif) = Ko{Di, D2), nor that if existing, it is unique. However, in the 
following important cases such an extremal mapping (for Kj or Kq) is known 
to exist: 

Theorem 6.6 (Gehring-Vaisala [21], Gehring [23]). The extremal mappings 
for Kj and Kq exist if 

(1) Di or D2 is a ball; 

(2) The boundary of Di, dDi has k components, where 2 <k < 00; 

(3) Di,D2 are tori in M^. 
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